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Abstract 

We explore the possibility of describing our universe with a singularity-free, 
closed, spatially homogeneous and isotropic cosmological model, using only general 
relativity and a suitable equation of state which produces an inflationary era. A 
phase transition to a radiation-dominated era occurs as a consequence of boundary 
conditions expressing the assumption that the temperature cannot exceed the 
Planck value. We find that over a broad range of initial conditions, the predicted 
value of the Hubble parameter is approximately 47 km- s _1 - Mpc -1 . Inflation is 
driven by a scalar field, which must be conformally coupled to the curvature if 
the Einstein equivalence principle has to be satisfied. The form of the scalar field 
potential is derived, instead of being assumed a priori. 



To appear in The Astrophysical Journal. 



1 



1 Introduction 



Singularity-free cosmological models have received renewed attention in the past few 
years (Rosen 1985, Israelit & Rosen 1989 hereafter IR, Balbinot et al. 1990, Starkovich 
& Cooperstock 1992 hereafter SC, Brandenberger et al. 1993). In the paper by IR, 
the universe was modelled as a closed Friedmann-Lemaitre-Robertson- Walker (FLRW) 
spacetime, bouncing from a Planck mass and radius state with inflation, emerging into 
the standard model of radiation and then matter dominance. The IR "prematter" 
heated with expansion as it emerged from the cold non-singular Planck state. This is an 
attractive simple picture, in contrast to the earlier prevailing theme among cosmologists 
that the universe emerged from a singular state in the big bang into a hot radiation- 
dominated universe, cooled to the point of vacuum energy dominance leading to the 
onset of inflation with further cooling. Then scenarios had to be constructed to re-heat 
the universe to the onset of the radiation-dominated phase of the standard model. 

SC built upon the IR model, incorporating a scalar field to describe all phases of 
the universe evolution as an integrated, albeit simplified, field theory. Besides being a 
physically meaningful source of gravitation in the early universe, a scalar field provides 
a perfect fluid stress-energy tensor also in the later epochs of the universe's history. 
Moreover, it has been proposed as a dark matter candidate (Starkovich 1992; Kofman 
et al. 1993; Delgado 1993; McDonald 1993). 

We adopt the Gliner (1966)-Markov (1982) picture of Planck limits to physical quan- 
tities and postulate that when the Planck temperature T p i is reached, the phase change 
from prematter (with equation of state P — (— 1 + j)p, where 7 is a small parameter) 
to radiation (P = p/3) occurs. 

In this paper, the work of SC is extended, refining the calculations of the evolution of 
the scale factor. We now determine the value of the present time as that corresponding to 
the evolution of the radiation temperature to 2.7 K. We then find that for a broad range 
of initial conditions (no "fine-tuning" requirement), the present value of the Hubble 
parameter is approximately 47 km- s _1 - Mpc -1 . 

The paper by SC is built upon the minimal coupling of the scalar field to the Ricci 
curvature. While minimal coupling is usually considered, it has been shown recently 
(Sonego & Faraoni 1993) that conformal coupling is the only acceptable possibility if 
the Einstein equivalence principle holds as applied to scalar wave propagation. Thus, 
it is of immediate interest to focus our cosmological theory on a conformally coupled 
scalar field. 

Instead of adopting a particular form of the scalar field potential (such as, e.g., 
A0 4 ) as is usually done, we start with the assumption that the equation of state is 
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inflationary, and derive the form of the potential. This approach reverses the procedure 
usually employed, but is perfectly legitimate, and has already been used in the literature 
(Lucchin & Matarrese 1985, Barrow 1990, Ellis & Madsen 1991, SC). Also, to be noted 
is a somewhat similar approach, deriving the scalar field potential from the observed 
spectrum of density perturbations (Hodges & Blumenthal 1990, Copeland et al. 1993, 
Turner 1993, Lidsey & Tavakol 1993). 

We regard as acceptable those cosmological models that merge satisfactorily with 
the standard model without singularities and without the need for fine-tuning of the 
parameters. We are able to impose the conditions which are deemed necessary and 
deduce the potential which achieves these aims. 

In sec. 2 we present the basic equations of the cosmological model with a scalar field. 
In sec. 3 we introduce boundary and initial conditions inspired by the Gliner-Markov 
ideas, and we solve analytically the equations for the dynamics of the universe. In sec. 4 
we solve numerically for the evolution of the scalar field and we find the form of the 
scalar field potential. The results are discussed in sec. 5. 



2 Basic equations 

Our theory is based on a closed, spatially homogeneous and isotropic FLRW model. The 
metric is given byQ 



ds 2 = dt 2 - a 2 (t) 



dr 



2 



+ r 2 (d6 2 + sin 2 9 d V 2 ) 



(2.1) 



where t is the comoving time, r G (0,1), 9 G (0, it), and tp G (0, 2ir). We assume 
that the only source of gravitation is a scalar field <j>(t) during the entire history of the 
universe. In the contemporary literature it is customary to assume that the dynamics 
of the universe are dominated by a scalar field during an early inflationary epoch. We 
extend this assumption to the later radiation- and matter-dominated epochs, in which 
a scalar field provides the simplest field-theoretic realization of a perfect fluid stress- 
energy tensor. However, the use of the scalar field as a source of gravitation in these late 
eras may be more meaningful; in fact it has been pointed out (Starkovich 1992, Kofman 
et al. 1993, Delgado 1993; McDonald 1993) that the 0-field could be a candidate for 
dark matter in the present era. 

1 We adopt the notations and conventions of Birrell & Davies (1982). 
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Most of the contemporary literature is built upon the minimal coupling of the scalar 
field to the Ricci curvature R. The Klein-Gordon equation with more general coupling 
is 

dV 

+ W + — = 0, (2.2) 

where 0<f> = g^V ^V^0, which assumes the value + 3<pa/a for our metric, £ is a 
numerical constant, and V{<fi) is the scalar field potential. While minimal coupling 
(£ = 0) is usually considered, other values of £ are possible. Recently, Sonego & Faraoni 
(1993) showed that the value £ = 1/6 ("conformal coupling") is the only acceptable 
possibility if scalar waves in a curved spacetime are to propagate locally as they do in 
Minkowski spacetime, as required by the Einstein equivalence principle. We stress that 
the derivation of this result is independent of the conformal structure of spacetime and 
of conformal transformations: in other words, conformal invariance of the Klein-Gordon 
equation is not required. In fact, eq. (|2.2j ) with any non-constant potential (or a mass 
term) is not conformally invariant^. 

Thus, we focus our cosmological theory on a conformally coupled scalar field. The 
stress-energy tensor is (see e.g. Birrell & Davies 1982, p. 87) 

= ^V^V^-^V^^V^^^-i^V^V^^-i^ 2 (r pv - -g tu ,R\+Vg fll/ +-<f>n<f>g ltt , . 

(2.3) 

This can be written in the form of a perfect fluid stress-energy tensor: 

T M „ = (p + P) u^u v - Pg^ , (2.4) 

where 

= 9 m (2.5) 

(V a 0V*0) 1/2 

is the four-velocity of comoving observers (we restrict ourselves to the consideration of 
a real scalar field which satisfies V Q 0V a > 0). The energy density and pressure are 
given by (a dot denotes differentiation with respect to t) 



P = T m = \ (0) 2 + \ 







a 




) 


+ - 






a 



+ V , (2.6) 



2 This can be understood physically by noting that the introduction of a mass, or a potential, sets a 
preferred length scale for the scalar field theory, which is then not invariant under a change of lengths 
and norm of vectors arising from a (position dependent) rescaling of the metric. Thus, conformal 
invariance follows from "conformal" coupling only in the absence of a potential. 
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i(p-T%) = i(0) 2 + i-00 + 
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3 a 



dV 1 
~d6 + 6 



7d N 


\ 2 11 











-v, 



(2.7) 



where eq. Q2.2|) has been used. It has been shown by Madsen (1988) that the stress- 
energy tensor of a scalar field can be put in the form corresponding to a perfect fluid 
only under the assumptions of spatial homogeneity and isotropy. If spatially dependent 
perturbations 6(j)(t, x) to the scalar field are considered, T^ u can still be put in the form 
corresponding to a fluid, but new terms appear, which can be interpreted as a heat 
flux and anisotropic stresses. These terms disappear in the case of minimal coupling 
(Madsen 1988). This has, as a consequence, the fact that, during a given era, there is no 
entropy production for the unperturbed field <f>(t). This is also true for spatially depen- 
dent perturbations of a minimally coupled scalar field, but, interestingly, it no longer 
holds for conformal coupling, where entropy production does occur. To our knowledge, 
perturbations of a non-minimally coupled scalar field have not been considered in the 
literature. We will not deal with this in the present paper, leaving it for future analysis. 
The Einstein equations 



-8ttT, 



give 



1 8tt 



a ( 2\ 

; + «*(7-3,)/> = 



where 7 is defined by the equation of state 

P=(l-l)p. 



(2.8) 

(2.9) 
(2.10) 

(2.11) 



For simplicity, we assume that 7 is constant during each era in which we divide the 
history of the universe. We will consider three eras, corresponding to different values 
of 7. We consider first the evolution of the scale factor, and then we will turn to the 
dynamics of the scalar field. 



3 The dynamics of the universe 

The evolution of the scale factor can be solved separately from the dynamics of the 
scalar field. The Einstein equations ( f2.9[ ) and ( |2.1(J ) give an equation involving only the 
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scale factor 

a /3 \ / a 2 + 1 



a + {2-'- l ){— J=°- 
In order to solve this equation, we introduce the conformal time rj defined by 

dt = a(r])dr] . (3.2) 

In the following, we will solve for the dynamics of the universe and of the scalar field in 
terms of the conformal time. The transformation to comoving time is given explicitly in 
eq. ( |3.43| ), and the final results are also given in comoving time. 
Equations fl210|) and Q become 

2 







fa' 








a" 




'a'\ 


a 3 


-( 


y) 


a" 

— + 

a 


(1 


7 - 



1 8tt 

+ ^ = T P, (3.3) 



' " I) P = , (3.4) 

where a prime denotes differentiation with respect to r\. If we focus our attention on the 
variable 

a' din a , 

u = - = — — , 3.6 
a arj 

we obtain, from eq. ( |3.5| ), the Riccati equation 

u + cu 2 + c = , (3.7) 

where 

c=\l~l. (3.8) 



In the following, we will consider values of 7 such that c 7^ 0. Equation ( p77|) is easily 
solved by setting 

1 w' 



u = 

c w 



ln(w 1/c )] , (3.9) 
which gives the real solution 

a(ri) = ao[cos(r]c + d)] 1 ^ , (3.10) 

where ao and d are integration constants. 

We consider three periods in the history of the universe: 
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• "prematter" (inflationary) era, < rj < i] r , in which 7 is small and positive 
(typical values are 7 P ~ 1CT 3 ). This gives an inflationary equation of state close 
to P = — p (the case 7 P = corresponding to eternal inflation is excluded). The 
smallness of the parameter 7 p determines the amount of inflation that the universe 
has experienced. 

• radiation era, r] r < r] < i] m , with 7 r = 4/3, corresponding to the equation of state 
P = p/3. 

• matter era, i] > i] m , with 7 m = 1, corresponding to a zero pressure dust 

(in the following, subscripts or superscripts p, r, m denote the prematter, radiation and 
matter era, respectively). In our model, the transitions between different eras, which 
occur at the conformal times i] r and i] m are imposed by discontinuous changes in the 
equations of state. The phase change at r\ r follows from the Gliner (1966)-Markov (1982) 
hypothesis of Planckian limits to physical quantities, in this case temperature, and the 
phase change at i] m occurs when the radiation density falls to the level of the non- 
relativistic matter density. Although the imposition of discontinuous phase changes is 
not derived from an analysis of the microphysics, the microphysics relating to the scalar 
field is itself unknown (Ellis 1991). However, the present procedure does yield the form of 
the scalar field potential V{<p), in contrast to the usual procedure of assuming, a priori, 
a form for V(<p). This is a variation of the approach first proposed by Synge (1955) and 
which has since been followed by various authors (Lucchin & Matarrese 1985, Barrow 
1990, Ellis & Madsen 1991, SC - see also Hodges & Blumenthal 1990 and Copeland et 
al. 1993). As applied to cosmology, this leads to the treatment of the early universe as 
a "hot laboratory for the particle physics relevant at these very early times" (Ellis & 
Madsen 1991). 

3.1 Boundary conditions and solutions for the scale factor 

Our initial and boundary conditions for the dynamics of the universe follow from the idea 
that the physical quantities are limited by the Planck values. This idea is present in the 
papers of Gliner (1966) and Markov (1982), and has been used by Rosen (1985), IR and 
SC. The closed cosmological model has a contracting phase preceding the turnaround, 
or "bounce" , when the Planck density is reached. We set rj = to be the time of the 
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bounce, at which point the scale factor has a vanishing derivative^, 

a'(0)=0, (3.11) 

p(0) = p pl , (3.12) 
(where p p i = 5.1566 ■ 10 93 g ■ cm -3 is the Planck density), which also imply 

a(0) = 4 p) = J— — , (3.13) 




p'(0) = . (3.14) 

At this stage, the universe is said to be in the "prematter" phase. Equation (|3.13| ) 
follows from eq. (|3TB|), and eq. Q3.14| ) follows from 



a 1 



p' + 3 7 -p = 0, (3.15) 
a 

which can be derived from the conservation equation VT^ = and the equation of 
state ( 2.11|) . The scale factor in the different eras is 

( [cos{ V c p )] 1/Cp 0<r]<r] r , 
a (v) = I at ] cos! 7 ? + d(r) ) Vr < V < r] m , (3.16) 

[ 4 m) cos 2 (r//2 + d {m ^) Vm<V ■ 

By imposing that the scale factor and its derivative (and hence, also the Hubble param- 
eter) are continuous at the transitions between the different eras (j] r and r) m ), we derive 
appropriate values for the integration constants: 

d {p) = , (3.17) 

aP = 0$ [(xnfocj,)]*- 1 , (3.18) 
d {r) = V r(cp-l) , (3.19) 

(r) 



a W = a ° - - , (3.20) 
cos(t^ + d«) ' v ; 

d (m) __ Vm _|_ d (r) _ ( 3 21 ) 



3 The possibility of imposing the initially static condition ( |3.11 ) is by no means trivial, and is guar- 
anteed by the inflationary equation of state that we assume, which circumvents the Hawking-Penrose 
singularity theorem. An equation of state satisfying the strong energy condition would imply a singular 
derivative of the scale factor, and an initial singularity. The inflationary equation of state violates the 
strong energy condition. 
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3.2 Other quantities 



We will now derive expressions for the conformal times r/ r and r/ m marking the duration 
of the prematter and radiation eras, respectively. In our model, as in IR and SC, the 
universe starts in a very cold state and heats during expansion in the inflationary era. 
This follows from the equation of state and, as a result, there is no requirement for 
"re-heating" as there is in other models of the early universe. We assume, according 
to the above mentioned idea of limiting values for the physical quantities, that inflation 
stops when the temperature reaches the Planck value T p i = 1.4169 • 10 32 K. The energy 
density at this time (r] r ) is 



p(Vr 



15 



(kT } 



pi j 



(3.22) 



where k = 1.38 -10 16 K 1 is the Boltzmann constant. 

Equations flp), fTI|) and (|3T22|) , together with (|3~T3|) , (fDJD and (pT9|) , give 



1 

rj r = — arccos 



15 p p i 



Cp/3y p 



7T Z 



(kr pl y 



(3.23) 



The dimensionless quantity in square brackets in eq. ( |3.23| ) will be used in the following, 
and its numerical value, with the velocity of light and the (reduced) Planck constant 
restored, is 

^' = 1.5201 . (3.24) 

We can derive an expression for r/ m by considering the evolution of the temperature T, 
which is described by the equation 



V a' 
7'— 7 — — 3 — 7 (7 — 1) = 



(3.25) 



T a 

(corresponding to eq. (2.20) in SC). The solution during the radiation era (7 = 4/3) is 

a(rj) T(r)) = const. = 0(77*) T(t}+) , (3.26) 

where 77* is any given value of 77 in the radiation era. In particular, for (77, 77*) = (r/ m , 77 r ) 
we get, from Eqs. ( |3.16| ) and (|3.19| ), 



arccos 



T, 



P i 



cos(rj r c p 



+ r/ r (l - c. 



VI ' 



(3.27) 
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where T m is the temperature at the last phase transition i] m . Kolb and Turner (1990, 
p. 77) give 

T m = T now • 2.32 • 10 4 n h 2 K , (3.28) 



where f2 = pj p c , Pc = 3H 2 / (8n) is the critical density, h = H / (100 Km • s 1 • Mpc 1 j , 
and T now = 2.7 K. We use the numerical value Q h 2 = 0.396 (which is compatible with 
a rather wide range of acceptable values for Q and h), giving T m = 2.4805 • 10 4 K. 
Actually, it will turn out that our final results depend only weakly on the value of T m , 
and are much more strongly affected by the choice of the parameter 7 P in the prematter 
era. 

We note that eq. ( |3.27|) sets a limit on the acceptable values of j p . Using eq. ( |3.23|) , 
the requirement that the argument of the arccos function in eq. (|3.27|) is not larger than 
1 gives 



1p< 3 



1 In {T p i/T„ 



2 In 1.5201 



(3.29) 



At this stage, it is useful to deduce the order of magnitude of the quantities involved, 
and make some approximations. We first note that the expression ( |3.16| ) for the scale 
factor in the prematter era becomes singular at i]\c p \ = tt/2. Typical values of 7 P that 
we consider are of order 10~ 3 , which gives 

AT = _ M ~ _300 (3.30) 

and, from eq. ( |3.23| ), 

7] r \c p \ = arccos(1.5201 JV ) = | + 8 = ~ - \S\ , (3.31) 

with \5\ ~ 10 -55 . Therefore, f] r \c p \ < tt/2, and inflation stops before the singularity in 
eq. (|3.16|) is reached, but numerically very close to tt/2. 



We now introduce 



x = ^\5\ (3.32) 



which, for T m ~ 10 4 K and 7 P ~ 10 3 , is of order 10 27 . These numbers justify the 
following expansions: We have 

cos(^ r c p ) = cos(7r/2 + 5) = — sin 5 ~ \5\ (3.33) 
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and eq. ( |3.27| ) gives 

r] m = arccos a; + r/ r (l — c p ) = —n/2 + x + r/ r (l — c p ) + 0(x 2 ) , (3.34) 

where we invert the function coss in the interval — ix < s < 0. Any other branch of 
arccos a; does not give the desired (increasing) behavior of the scale factor after the last 
phase transition. 

We now compute the age of the universe, f] now , in conformal time. To this end, 
we note that after the transition between the radiation and matter eras, radiation is 
decoupled from matter, and behaves like a perfect fluid with 7 = 4/3. Equation ( |3.26| ) 
gives, for (77,77*) = (r]m,Vnow) 

0(77771) T m = a{jinovu) T now (3.35) 

Vnow j(m)\ _ 2 ( Vn 



and using eq. (|3.16| ) 

2 



cos 



Now, Eqs. (|3~3l) , (|3~T9| ) and (|3~2l|) give 



n 

and therefore 



-^ + d (m) -arccosa; (3.37) 



Vnow + d(m) \ = T pl ^ (3 ^ 3g) 



COS , . f -^— — 

where we have chosen the positive sign. We must now invert the function coss in the 
interval (— it, 0) as before. This branch (which we call "branch 2") is the opposite of 
that obtained by inverting coss in the interval (0,7r) (denoted by "branch 1"). Any 
branch other than branch 2 will not give the desired behavior of the scale factor (i.e. 
increasing for r] m < r] < r] now , and for r] now < r) < rj 

maxi where f] max is the time of 
maximal expansion of the universe before it starts recollapsing) . Therefore, we write 



f]now — 2 arccos 



± m \ 

■ X 



2.7K 



-rim- 2d^ . (3.39) 

branch 1 



Equation ( 3.39Q gives an upper limit on j p . The requirement that the argument of the 



arccos function in eq. ( |3.39p is not larger than 1 gives, using eqs. ( |3.30| )-( p.33| ) 



In 1.5201 

1p < 7 , . 3.40 

3 In [Vl^201T pl /y/T m T nc ' 
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The conformal time of maximum expansion of the universe is trivially given by eq. ( |3.16| ) 

as 

Vmax = -2d™ ; (3.41) 



it is reached when 



a(r] max ) = a { ™> . (3.42) 

It is useful to translate conformal times into comoving times. Conformal time is 
defined by eq. (|3.2|) . Choosing t = at r\ = we get, from eq. ( ft.lOQ , 



t{rf) = ao dr\ [cos(r/c + d)] 



l/c 



(3.43) 



This integral has to be computed numerically in the prematter era, while in the radiation 
(c r = 1) and matter (c m = 1/2) eras, the integration is trivial. The expressions for the 
comoving times corresponding to rj r , r) m , r\ now and r] max are: 



» 



Vr 



drj 



1 



t r + a. 



(r) 



o [cos(r/c p )] 1/|Cpl ' 
sin(r? m + d {r) ) - sm(r) r c p ) 



(m) 



v now,max "m 



tin ~l~ 



Vnow,max Vm ^^-{V now, max 2t£ ) sin(^ m 4" 2(^ 

The Hubble parameter at the present time is 

a'(r] now ) 9.2503 ■ 10 29 sm(r) now /2 + d^ 



m) "> 



H(r] n 



Km ■ s^ 1 ■ Mpc 



-i 



(3.44) 

(3.45) 
(3.46) 

(3.47) 



~a 2 (w,)" a { m) cos 3 (r/ nou ,/2 + rfM) 

The evolution of the energy density is governed by eq. ( |3.15| ), which has the solution 

pa 37 = constant = p(r/*) a 37 (77*) , (3.48) 

where 7]* is any fixed conformal instant during a given era. By choosing (?7*, 77) = (0, rj r ), 
mi T]now ), {r]m,Vmax) we get, respectively, 



P(Vr 



P P i 



a(r] r 



37p 



p(Vm) = p{Vr) 



airjr 

a{Vrr 



1 4 



(3.49) 
(3.50) 
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PiVnow) = PiVm) 
p(Vmax) = P(j] m ) 



a(rj r , 



(m) 



(3.51) 
(3.52) 



3.3 Numerical values 



We present tables with numerical values for the conformal and comoving times marking 
the transitions between the prematter, radiation and matter eras, for the present age of 
the universe and for the time of its maximal expansion when it just starts recollapsing. 
We also give the corresponding values for the scale factor (in cm) and energy density (in g- 
cm -3 ). The present value of the Hubble parameter (in Km-s _1 -Mpc _1 ) is also computed. 
The tables correspond to the following choices of parameters: T m = 2.4805 • 10 4 K and 
7 P = 1.8500 • 10" 3 , 1.9000 • 10~ 3 , 1.9500 • 10" 3 , 2.0000 • 10" 3 , 2.0100 • 10~ 3 , 2.0153 • 10" 3 , 
which are allowed by the limits in eqs. ( |3.29p and fl3.40| ). It is to be noted that for a 
broad range of initial conditions corresponding to values of 7 P between 1.75 • 10~ 3 and 
1.95 • 10~ 3 , the value of the Hubble parameter hardly varies, remaining close to 46.7 
km- s _1 - Mpc -1 . If the preferred model of the universe is one which does not demand 
the fine-tuning of parameters, it could be argued that this is the proper value of H . 
While this value is relatively low, it is compatible with the currently accepted range. 

It is to be noted that in the narrow range of allowed values of 7 P near the critical point 
from 2.01 • 10 -3 to 2.0323 ■ 10~ 3 , the predicted H drops rapidly beyond the acceptable 
range. Moreover, while the value of H is substantively governed by the choice of the 
parameter j p , it depends only weakly on the value of the parameter T m . 



4 The scalar field dynamics 

We turn now to the dynamics of the conformally coupled scalar field driving inflation. 
As discussed above, we do not assume a priori a form for the scalar field potential, but 
rather we derive it as a consequence of the assumption that the equation of state changes 
in the different eras of the universe. The Klein-Gordon equation ( [2.2|) with £ = 1/6 is 

4 > » + 2-^+(- + l) ( p + a 2 ^- = 0. (4.1) 
a \ a J acp 
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Equations QSTTlD , ( pTFQ and Q2^) give 



P 2 . ., 2 4 a' 

7 = - + 1= 2 - 0' +--• 
p 3 3 a 



1 2 dV 2 , 2 



' + 2a 2 V + 




(4.2) 



from which we get 

v 1 /4 




2 7 a 2 V3 

Now, using eqs. (|2.6|) and (|3.2|) in eq. Q2.9Q we find 



i + 



— <j) — 
37 rf^ 



(4.3) 



Substituting the value of dV/d(f) given by eq. ( |4.1| ) into eq. 
eq. flOl), we find 



(4.4) 



and the result into 



2l. 
a 



2(0') 2 + 2 





97 


[(+ 


- 1 


a \ a / 


+ 87 







. 



(4.5) 



The last equation can be solved numerically in the prematter era, once the parameter 
7 P is fixed, and initial conditions <po, <fr' are specified. Then, one can solve numerically 
for the values of <fi'{rj) and V(rj) from eq. ( |4.4| ) and, by eliminating the conformal time, 
one obtains a numerical form for the scalar field potential V(<f>). Using the MAPLE 
computer program, we have done this for 7 P = 2.0153 ■ 10~ 3 (corresponding to Table 6). 
The result is plotted in two different ranges of values of (p. The resulting V(<p) does not 
depend on the particular solution <f>(rf) coming from specific initial data (</> , <fi' ), as must 
be the case. Integration gives only narrow ranges of (j>{rj) for given initial conditions, so 
the relation - V is plotted in narrow windows spanning only a part of the possible 
range of values of (p. We treat the initial conditions as arbitrary. The question if there 
are preferred initial conditions, and then a corresponding particular range for <fr remains 
an open problem. 

The features of the resulting V(<f)) are as follows: the scalar field potential has the 
symmetry V(<j)) = V(— <p), as to be expected from an examination of eq. (|4.4|), and for 
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4> > the overall shape is plotted in fig. 1. V is slightly positive for small values of 
0, and becomes negative as increases. Superposed upon the basic overall shape are 
variations of V(<f>) on a small scale, which occur in narrow ranges of 0; some of these 
variations, with local maxima and minima are plotted in fig. 2. 

If one accepts the above mentioned idea that the scalar field is the dominant com- 
ponent of dark matter also in the post-inflationary eras of the universe, it is immediate 
to derive the form of the potential V in the radiation era. In fact, the trace of the 
energy-momentum tensor is 

T a a = p-3P = W-<P^, (4.6) 
which vanishes when 7 = 4/3, giving 

V(4>) = A0 4 , (4.7) 

where A is a constant. This form of the scalar field potential has been used widely in 
the literature. It should be noted that this result holds only in the case of conformal 
coupling. For example, when £ = 0, SC obtained a different form of V. 



5 Discussion and conclusion 

We have constructed a closed cosmological model which has no singularities. In contrast 
to other works (e. g. Balbinot et al. 1990, Brandenberger et al. 1993), we use only 
general relativity, and an inflationary equation of state in the early universe. We impose 
discontinuous changes in the equation of state, and solve for the dynamics of the universe. 
The key ingredient to obtain a singularity-free model is the Gliner-Markov idea that the 
physical quantities are limited by the Planck values. Our boundary and initial conditions 
on the scale factor reflect this idea. The evolution equations for the scale factor and 
the energy density can be solved exactly by using conformal time instead of the more 
commonly used comoving time. However, the physically significant conformal times have 
been converted to comoving times. We considered three different eras for the history 
of the universe, corresponding to different values of the constant 7 in the equation of 
state. The early inflationary era is dominated by prematter, then a radiation-dominated 
and a matter-dominated era follow. We determined the times of transition between the 
different eras, the present age of the universe, the time of its maximal expansion, and 
the Hubble parameter, for various values of the parameter j p . We obtained values of 
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H and t now in a range which agrees with the observations and the current theoretical 
expectations. 

The only source of gravitation is provided by a scalar field, which is chosen to be 
conformally coupled to the Ricci curvature, for the above mentioned reasons. We pro- 
ceeded to the numerical solution of the equations for the scalar field and the potential as 
functions of conformal time. From these, we derived the form of the potential V((f>). It 
is interesting to note that the derived potential differs from that computed in SC for a 
minimally coupled field. Moreover, the greater complexity of the non-minimal coupling 
prevents one from deriving a differential equation, and then an analytical expression for 
V(4>), as was possible in the SC paper where the coupling was minimal. 

The potential V(4>) exhibits variations on a small scale, with local maxima and 
minima. One could ask which range of values of <fi is relevant to the physics of the early 
universe, and hence what is the shape of the potential for the appropriate range of the 
scalar field. The answer to this problem is not an easy one, since the actual values 
of (ft are determined by the arbitrary initial data set (<f)o,<f)' ). Although one can try 
to constrain the value of <fro with arguments based on mass scales (a very speculative 
approach, but one which has received some attention in the literature - see e.g. Hosotani 
1985 and Futamase & Maeda 1989), there are no indications on how to constrain the 
initial derivative (fi' . It turns out that the minimally coupled case requires only the 
consideration of one initial value, 0o- 

It must be noted that various inflationary scenarios with a nonminimally coupled 
scalar field and an a priori chosen form of the potential encounter problems (Abbott 
1981; Starobinsky 1981; Futamase & Maeda 1989; Futamase et al. 1989; Bruni 1993). 
On the other hand, conformal coupling is required if the Einstein equivalence principle is 
to hold. In our approach, because of the fact that we do not fix a form of the potential, 
but rather derive it, the previous difficulties are avoided, and we can have both conformal 
coupling and a viable inflation (at least to the extent explored in this paper). 

Our model universe starts very cold and initially static from Planck size and density, 
and undergoes an inflationary expansion, during which it heats. It must be noted that 
the possibility of imposing the initial condition a'(0) = 0, and the absence of an initial 
singularity are by no means trivial, in the light of the Hawking-Penrose singularity 
theorems (see the formulations in Hawking & Ellis 1973; Wald 1984; and references 
therein). In fact, the strong energy condition required by these theorems, which can be 
written as 

p + P>0 and p + 3P>0, (5.1) 
is violated for 7 < 2/3, which is our case, since 7 P ~ 10~ 3 . On the contrary, in the 
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standard model starting from the big bang in a radiation-dominated epoch, 7 = 4/3. 
The weak (p > and p + P > 0) and dominant (p > \P\) energy conditions are satisfied. 
The inflationary equation of state permits the avoidance of singularities, and in general 
relativity, inflation (defined as positive acceleration of the scale factor) occurs if and only 
if the condition p + 3P > is violated. 

Regarding the regularity of the functions of physical interest, the continuity of the 
scale factor a(rj) and its derivative were assured by the proper choice of integration 
constants in eqs. (|3.17|) — (|3.21|) . Thus, the Lichnerowicz (1955) conditions on the metric 
are satisfied and from eq. ( |2.9|) , the density p is continuous at the phase transitions. 
Since we choose different constant values of 7 at the three different universe epochs 
(j p ~ 10~ 3 , 7 r = 4/3, 7 m = 1 respectively), the pressure from eq. (|2.11|) is discontinuous 
at the phase transitions. 

Since the idea of treating the scalar field as the primary source of gravitation in the 
radiation- and matter-dominated epochs is a speculation at this stage, we only found 
V{4>) in the radiation era (eq. ( |4.7| )) and we did not present a numerical solution for <fi 
and V{4>) otherwise for those epochs. We have treated the initial data (</>o,</>o) f° r the 
scalar field as arbitrary and we can use the values of <fi and <$' at the end of a given epoch 
as the initial values for the succeeding epoch. Thus, it is always possible to evolve a 
continuous solution in 0, 0'. However, since the potential V(<p) has a different functional 
form in the different epochs, corresponding to the different equations of state, it is not 
continuous. The difference in form is evident in comparing the numerical solution in 
figs. 1, 2 with the analytical form in eq. ( f4.7|) . However, if 0(r?), 4>'{t]) are continuous, 
then the potential as a function of the conformal time, V(t]), will be continuous as can 
be seen in eq. ( |2.6|) . However, this is not inconsistent with the discontinuity of V as 
expressed function of d>. 

After the universe reaches the point of maximal expansion, it experiences its past 
history again in reversed order, until it reaches a point, in a prematter era, in which 
the derivative of the scale factor vanishes, and then the past history is repeated. The 
vanishing of a' makes it possible to match the two FLRW metrics before and after 
the bounce in such a way that the scale factor and its derivative are continuous, thus 
satisfying the Lichnerowicz (1955) conditions. 

We did not address the problem of density perturbations and the formation of struc- 
tures in our model universe. Perturbations of a non-minimally coupled scalar field suffer 
from the problems pointed out in sec. 2, and will be the subject of future analysis. Possi- 
ble generalizations of this paper include the study of open models and of inhomogeneous 
metrics. 
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Figure 1: the overall shape of the scalar field potential V((f>) (rescaled by 10 68 cm -2 , 
in units such that G = c = 1) for > 0. The plot is obtained by superposing different 
windows like that in fig. 2, each spanning a small range of values of <fi. The plots in fig. 1 
and Fig. 2 correspond to T m = 2.4805 • 10 4 K and 7 P = 2.0153 • 10~ 3 . 

Figure 2: an example of the variations of V(4>) on a small scale, with local maxima 
and minima. 
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